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METHOD AND ARRANGEMENTS REIiATING TO SATELLITE - BASED POSITIONING 
TECHNICAL FIELD OF THE INVENTION 

The present invention relates in general to positioning of 
mobile equipment by use of satellites and in particular to such 
positioning assisted by land based communication nodes. 

DESCRIPTION OF RELATED ART 

In recent years, determination of the geographic position of an 
object, equipment or a person carrying tlie equipment has become 
more and more interesting in many fields of application. One 
approach to solve the positioning is to use signals emitted from 
satellites to determine a position. Well-known examples of such 
systems are the Global Positioning System (GPS) (see e.g. [1]) 
and the coming GALILEO system. The position is given with 
respect to a specified coordinate system as a 
triangulation/trilateration based on a plurality of received 
satellite signals. 

A stand-alone GPS receiver can obtain, full locking to GPS 
satellite signals, without having any other information about 
the system except nominal carrier frequency and the rules by 
which data carried by the signals are modulated. Basically, the 
three-dimensional position as well as a receiver clock bias to 
the satellite time have to be determined in the position 
calculation step. 

Assisted GPS (AGPS) has been defined as an enhancement of GPS 
(see e.g. 3rd Generation Partner- ship Project (3GPP) 
specifications TS 25.331 or TS 44.031 or Open Mobile Alliance 
(OMA) specifications for Secure User Pla.xie Location (SUPL) ) for 
integration of GPS receivers into user equipment, i.e. mobile 
stations, of cellular communication systems. Assisted GPS in 
general aims at improving the performance of GPS receivers in 
many different respects, including detection sensitivity, time 
to obtain a location estimate, accuracy and saving battery 
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power. This is done by moving some functionality from the GPS 
receiver in the mobile station to the network and hence only 
performing a subset of the GPS tasks in the GPS receiver itself. 

There are two types of AGPS, Mobile Station (or User Equipment) 
5 based and Mobile Station (or User Equipment) assisted. In Mobile 
Station based AGPS, the location of a mobile station is 
calculated in the mobile station using ranging signal 
measurement results determined by the mobile station together 
with assistance data provided by the network. In Mobile Station 

10 assisted AGPS (sometimes also referred to as Network based 
AGPS) , the mobile station only measures and reports timing of 
received ranging signals reflecting the pseudoranges to the 
Space Vehicles (i.e. satellites). For both types of AGPS, the 
measured timing of the ranging signals are truncated modulo 1 ms 

15 which corresponds to a distance of 3 00km. When calculating the 
mobile station location, either in the mobile station itself or 
in a network location server, the complete pseudoranges need to 
be reconstructed using apriori information about the mobile 
station location together with the ranging signal measurement 

20 results determined by the mobile station, in order to compute 
the precise mobile station location. 

The inventor of the present invention have identified a problem 
with AGPS in that measuring and reporting truncated timing of 
received ranging signals may cause ambiguities when determining 

25 pseudoranges to space vehicles, if the precision of the apriori 
information about the mobile station position is too low, i.e. 
the uncertainty of the mobile station position in said apriori 
information is too large. As a result, if an incorrect 
pseudorange is selected and used as a basis for determining the 

30 location of the mobile station, there will be a significant 
error in the calculated mobile station position in the order of 
e.g. 100 km. 
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SUMMARY OF THE INVENTION 

The problem dealt with by the present invention is providing 
increased robustness, in the context of satellite based 
positioning with assistance data, against ambiguous pseudorange 
5 reconstruction. 

The problem is solved by a method according to claim 1, an 
apparatus according to claim 19 and a computer program embodied 
on a computer-readable medium according to claim 33. 

An advantage afforded by the invention is increased robustness 
10 against ambiguous pseudorange reconstruction in connection with 
satellite based positioning with assistance data such as 
Assisted GPS (AGPS) . 

Another advantage of the invention is that the increased 
robustness is achieved without reducing the detection 
15 sensitivity. 

Yet another advantage of the invention is that the increased 
robustness is achieved with an insignificant increase in 
processing delays * 

The invention will now be described in more detail with 
2 0 reference to exemplary embodiments thereof and also with 
reference to the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig 1 is a schematic view of an example scenario of Mobile 
Station assisted AGPS in which the present invention is applied 



25 



Fig. 2 is a diagram illustrating C/A code and navigation data 
bits in GPS ranging signals 

Fig. 3 is a block diagram illustrating the format of GPS 
navigation data in GPS ranging signals 
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Fig. 4 is a diagram illustrating time in different parts of the 
system illustrated in Fig. 1 

Fig. 5 is a flow diagram illustrating a basic method according 
to the invention. 

5 Fig. 6 is a schematic block diagram of a mobile station. 

Fig. 7 is a flow diagram illustrating processing performed by 
the mobile station of Fig. 6. 

Fig. 8 is a schematic block diagram of a location server 
providing an example embodiment of an apparatus according to the 
10 invention. 

Fig. 9A-B is a flow diagram illustrating a detailed example 
embodiment of a method according to the invention 

Fig. 10 is a diagram illustrating results and intermediate 
results when performing the method of Fig. 9A-9B. 

15 Fig. 11 shows an example of a computer-readable medium 

DETAILED DESCRIPTION OF THE EMBODIMENTS 

Fig. 1 illustrates a non-limiting example scenario in which the 
present invention may be applied. In this example scenario a 
basic wireless communication system SYSl together with the 
2 0 Global Positioning System (GPS) is used to provide Mobile 
Station assisted AGPS. The example wireless communication system 
SYSl illustrated in Fig. 1 is a Universal Mobile 
Telecommunication System (UMTS) . The communication system SYSl 
includes a network part NET1 and User Equipment (UE) , 

2 5 alternatively referred to as mobile stations (MS) . The network 

part NET1 comprises a core network CNl and a UMTS Terrestrial 
Radio Access Network (UTRAN) RANI. The core network CNl includes 
a Mobile services Switching Center (MSC) node MSC1 that provides 
circuit-switched services and a General Packet Radio Service 

3 0 (GPRS) node SGSNl , sometimes referred to as a Serving GPRS 
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Support node (SGSN) , which is tailored to provide packet- 
switched type services. 

Each of the core network nodes MSC1 and SGSNl connects to the 
the radio access network RANI over a radio access network 
5 interface referred to as the Iu interface. The radio access 
network RANI includes one or more radio network controllers 
(RNCs) . For sake of simplicity, the radio access network RANI of 
Fig. 1 is shown with only one radio network controller node 
RNC1 . Each radio network controller is connected to and controls 

10 a plurality of radio base stations (RBSs) . For example, and 
again for sake of simplicity, Fig. 1 only illustrates a first 
radio base station node RBS1 and a second radio base station 
node RBS2 connected to the radio network controller node RNC1 . 
The interface between the radio network controller RNC1 and the 

15 base stations RBS1 and RBS2 is referred to as the Iub interface. 
Mobile stations, such as mobile station MSI shown in Fig. 1, 
communicate with one or more radio base stations RBS1-RBS2 over 
a radio or air interface referred to as the Uu interface. 
Each of the radio interface Uu, the lu interface and the Iub 

2 0 interface are shown by dashed lines in Fig. 1. 

In Fig. 1, the GPS system is represented by Space Vehicles, i.e. 
satellites, SV1-SV4. Each Space Vehicle SV1-SV4 transmits a 
corresponding ranging signal RS1-RS4 . Please note that for sake 
of simplicity, only four Space Vehicles SV1-SV4 are illustrated 

2 5 in Fig. 1. 

When determining the position of the mobile station MSI in Fig. 
1 using mobile station assisted AGPS, the mobile station MSI 
receives assistance data from and reports measurement results to 
a location server 101. Based on the reported measurement results 

3 0 and apriori information on where the mobile station is located, 

the location server calculates the location of the mobile 
station MSI. Depending on how a location server is connected to 
a cellular network, AGPS can be divided into two categories, 
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"control plane solutions to AGPS" and "user plane solutions to 
AGPS" . 

In "control plane solutions to AGPS", the location server 
functionality (which may be implemented in a separate location 
5 server node, sometimes referred to as a Serving Mobile Location 
Center (SMLC) or Standalone SMLC (SAS) , or integrated together 
with other functionality in other network nodes such as radio 
network controllers) is tightly integrated with the cellular 
network and assistance data and measurement results are 

10 communicated using so called control plane signaling. This 
solution is further characterized in that typically the location 
server would receive information of in which cell a mobile 
station is currently operating and the location server would 
^PPly this information as the apriori location of the mobile 

15 station when calculating the location of the mobile station. 
Hence the uncertainty in the apriori location information 
corresponds to the cell size. 

In "user plane solutions to AGPS", the location server 
functionality is less closely integrated with the cellular 

2 0 network and assistance data and measurement results are 

communicated using so called user plane signaling, i.e. ordinary 
user data packets are used to convey this information 
transparently to the cellular network. This solution is further 
characterized in that the location server may not receive 
25 information of in which cell a mobile station is located or at 
least may not always be able to associate a given cell identity 
with a specific geographical area corresponding to the area 
covered by the cell. Hence, for user plane solutions to AGPS, 
the uncertainty in the mobile station apriori location 

3 0 information may be significantly larger than the cell size and 

may correspond to e.g. the size of the country in which the 
mobile station is currently operating. 



In the example scenario of Fig. 1, a user plane solution to AGPS 
is illustrated where the location server 101 is connected to the 
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cellular network NETl via an Internet Protocol (IP) based packet 
data network 102 . 

The GPS Space Vehicles SV1-SV4 transmit ranging signals RS1-RS4 
with a spectrum centered at 1575.42 MHz. Fig. 2 illustrates how 
each ranging signal RS1-RS4 includes a stream 201 of - navigation 
data bits 202 that are spread by a spreading code defined by a 
so-called Coarse/Acquisition (C/A) code 203 that is unique for 
the Space Vehicle transmitting the signal. The C/A code 203 has 
a length of 1023 chips and a chip duration of 1/I.023xl0 6 s, 
i.e. the C/A code comprises a sequence of +/- 1 that changes at 
a rate of 1.023xl0 6 Hz and repeats itself every 1 ms . The 
navigation bits 202 have a bit period of 20 ms, i.e. 
corresponding to 20 C/A code repetitions. 

The navigation data includes among other things a set of so- 
called ephemeris parameters that enables the receiver to 
calculate the precise position of the satellites at the time of 
signal transmission. The precise time of transmission can also 
be read from the navigation data. 

Fig. 3 illustrates more in detail how the navigation data is 
further divided into 5 subframes 3 01-305, each of length 6 
seconds. Each subframe 301-305 is divided into 10 words each of 
length 0.6 seconds and containing 3 0 data bits. A time stamp, 
GPS Time Of Week (TOW) , is transmitted in the second word, the 
Handover Word (HOW), of every subframe 301-305. The indicated 
time is the time of transmission at the end of the subframe in 
question. The TOW is thus repeated every 6 seconds. 

Each ranging signal RS1-RS4 basically defines a clock which is 
measured by the mobile station MSI. The clock indicates the time 
of signal transmission. If the mobile station MSI knows the GPS 
system time, then the clock reading can directly be used bo 
determine the time delay , and hence the range from the Space 
Vehicle transmitting the ranging signal to the mobile station 
MSI . By measuring three ranges and utilizing the knowledge about 
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Space Vehicle locations at time of transmission, the location of 
the mobile station MSI in three dimension can then be 
determined. However normally the mobile station MSI does not 
have knowledge about precise GPS system time, so one more 
5 measurement is needed to eliminate the mobile station clock 
bias. 

The sequences of Fig. 4 illustrate clock relations (expressed in 
milliseconds) for different parts of the system illustrated in 
Fig. 1. Each Space Vehicle SV1-SV4 carry precise atomic clocks 

10 to maintain clock stability. The Space Vehicle transmissions are 
however not perfectly synchronized to GPS system time as 
illustrated in Fig. 4. In Fig. 4, sequence 401 represents the 
GPS system time, sequence 411 represents the clock of Space 
Vehicle 1, sequence 41N represents the clock of Space Vehicle N, 

15 sequence 402 represents the clock of mobile station MSI of Fig. 
1 while sequences 421 and 42N respectively represent time as 
read in the ranging signals from Space Vehicle 1 and Space 
Vehicle N respectively received by the mobile station MSI. By 
drawing a vertical line 431 through the timing diagram one may 

2 0 obtain a snapshot of all clock readings as observed in various 
points in space. GPS system time 401 is defined as an ensemble 
average based on a set of ground station clocks and a subset of 
Space Vehicle clocks. As demonstrated in Fig. 4,- the individual ^ 
Space Vehicle clocks 411 and 41N and the mobile station clock 

25 402 are slightly offset (see SV clock biases 412 and 413 and 

mobile station clock bias 414 respectively) compared to GPS ^ 
system time 401. A model for the individual offsets of the Space 
Vehicle clocks is transmitted as part of the navigation message 
from each Space Vehicle. When the signals reach a point on the 

30 earth surface (e.g. the current location of mobile station MSI), 
they have been delayed with an amount depending on the range 
from the Space Vehicle in question to said point, on the earth 
surface. The delay is typically 60-85 milliseconds (ms) as 
illustrated by the clock readings in Fig. 4. 
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When determining the position of a mobile station using AGPS, 
the mobile station measures the C/A code boundary locations 
position in time with respect to a selected point in time, i.e. 
the C/A code phase, for the received ranging signals. The C/A 
5 code phase is determined modulo 1 ms (i.e. one C/A code period) . 

A mobile station implementing mobile station based AGPS 
calculates its position at the selected point in time based on 
the measured C/A code phases of the received ranging signals and 
assistance data (received from the network) including space 
10 vehicle ephemeris and clock correction data together with 
apriori information about the mobile station location. 

A mobile station, such as mobile station MSI in Fig. 1, 
implementing mobile station assisted AGPS instead transmits a 
wireless signal reporting the C/A code phases (expressed in 

15 terms of whole and fractional chips of the C/A code from the 
selected point in time until the beginning of the next C/A code 
repetition) for the received ranging signals together with an 
estimate of the GPS system time corresponding to the selected 
point in time. Based on the inf ormation reported by the mobile 

2 0 station and apriori information about the mobile station 
location, a location server, such as location server 101 in Fig. 
1, either in the cellular network or in another network 
calculates the position of the mobile station. 

The inventor of the present invention has recognized that the 
25 AGPS way of measuring the C/A code phases modulo 1 ms and hence 
characterizing each ranging signal by time mod 1 ms, causes 
problems when the uncertainty in the apriori information of the 
mobile station location is too large. As demonstrated below, 
when the initial location uncertainty is more than 75 km, the so 
30 called pseudorange to a space vehicle can not be unambiguously 
reconstructed. 

A GPS receiver basically measures the pseudorange to a number of 
satellites (please note however that in mobile assisted GPS, the 
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GPS receiver integrated in the mobile station does not perform 
the complete calculations needed to determine the pseudoranges, 
but only provides basic data needed to calculate pseudoranges) . 
The pseudorange is 

5 Pi = c- (t u -t ti ) t (1) 

where t u is the GPS receiver (integrated in the mobile station) 
clock reading at the time of reception, and t t i is the time of 
signal transmission of the ith Space Vehicle, and c is the wave 
propagation speed. The pseudorange differs from true range with 

10 a number of perturbing factors (receiver clock bias, ionospheric 
and tropospheric delays, Space Vehicle clock bias, measurement 
errors, etc) . For purpose of clarity, the influence of most of 
these error sources are neglected in the presentation below. 
There are known techniques to compensate for many of the above 

15 listed error sources (see e.g. [1], [2]). Furthermore the effect 
of Space Vehicle movement and earth rotation is also omitted as 
it is also well known in the art how to handle these effects 
(see e.g. [1], [2]). The simplified model is then that the 
measured pseuodrange obeys 

20 



Pi = I x u -x si J +b+Gi (2) 

Here x u = (x u y u z u ) is a row vector containing the three- 
dimensional coordinates of the unknown receiver location. 
Similarly x S i is the row vector containing the coordinates of 

25 the ith Space Vehicle. The notation \z\ means the norm of the 
vector quantity within brackets, which is equal to {zz V /2 . In 
this case it can be interpreted as the distance between the GPS 
receiver /mobile station and the Space Vehicle. Furthermore b is 
the receiver clock bias (expressed as a range) , 

30 b=c- (t u -t GPS ) (3) 
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where t GPS stands for GPS system time. Finally e± is the 
measurement error. 

The mobile station in mobile station assisted AGPS reports only 
pseudoranges modulo 1 C/A code period, ie the integer number of 
code periods in the pseudorange (1) is not known. It has to be 
reconstructed. Let R in general represent the range 
corresponding to the truncation interval used when measuring 
ranging signal timing. As an example, the range R corresponding 
to a truncation interval of one C/A code period (1 ms) is then 

R = c-lCT 3 (4) 

We thus have that 

Pi = k±R+v ± (5) 

where v± is the reported measurement that satisfies 0< v± < R, 
and where admissible values of the integer k± needs to be 
reconstructed. When the apriori initial location is large, 
several values of k i may be admissible. Furthermore the presence 
of a common bias term b makes it difficult to reconstruct k. 
exactly. However, from position calculation point of view it is 
the relative pseudoranges that matter, since any constant bias 
term will cancel in the position calculation. 

Define the reconstructed pseudorange as 

Pi*= ks+R+Vi (6) 

Next it will be shown that it is not necessary that k±*=k± for 
all i . However the relative pseudoranges must be correctly 
reconstructed, i.e. 

k±*=k± + X, for all i, X an integer (7) 

This can be explained as follows. Inserting (7) in (6) and using 
(5) , (2) gives that 
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pi*= k±R+XR+Vi =pi+XR =lx u -x sl l+b+XR+e± = jx u -x si j +b*+ei (8) 

i.e. that the reconstructed measurement equation (8) has the 
same structure as (2), the only difference being that b is 
replaced by b*=b+XR. 

5 The reconstruction is done in the following way. It is assumed 
that an apriori location x u0 of the mobile station is known, 
along with an uncertainty A, such that 

I x u -x u0 1 <A (9) 

First we determine the predicted pseudorange to Space Vehicle #1 
10 (arbitrarily selected among the Space Vehicles from which 
ranging signals have been received/measured) 

Pi ' = lx sl ~x u0 l (10) 

and next we try to find Jci* that satisfies 

k 1 *R+v 1 =p 1 ' (ID 

15 In most cases no exact solution is available so we round the 
result to the nearest integer so that 

Jcj* = round ( (pi'- Vi) /R) (12) 

So we get 

p x *= round ( (p! 1 - v^/RJR+Vi (12b) 

20 Next step is to estimate the set of admissible values k±* , for 
i>2. Admissible pseudoranges p±*, satisfy 

Pi*~Pl *=pi-Pl= /Xu-Xsi I ~ /Xu-Xsl I ( 13 ) 

In (13), measurement errors have been neglected. To obtain a 
sharp bound one may calculate the expected range difference 
25 (i\ e. difference between ranges to Space Vehicle #1 and Space 
Vehicle #i respectively) at all possible receiver locations that 
satisfy (9) , i.e. at all possible receiver locations within the 
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initial location uncertainty area defined by x u o and A. This way 
we obtain an inequality 

5p ± -A ± < p ± *-pi* < Spi+Ai , l=2,„..,n (14) 

where 8p± is the mean value of the maximum and minimum range 
5 difference that can be found anywhere within the initial 
uncertainty area. Similarly A± is half the difference between 
the maximum and minimum range difference that can be found 
anywhere within the initial uncertainty area. 

All admissible pseudoranges p±* , which in combination with pi * 
10 form admissible relative pseudoranges p±*-pi* , satisfy equation 
(14) . 

Using (6) , equation (14) can be rewritten as 

8p± -A ± ^ki+R+Vi-p!* <5pi+Ai (15) 

Admissible values for p±* can now be found by checking for 
15 varying ki*, whether the inequality (15) is satisfied. Note 
that several possible values of ki * may be admissible. Possible 
values of the reconstructed pseudoranges (6) are stepped with a 
spacing of R meters. The width of the uncertainty region is 
2 A± . So whenever 

20 2Ai > R (16) 

there is a risk that the reconstruction cannot be done 
unambiguously . 

A few alternative approaches for determining the admissible 
pseudoranges is presented next. 

25 It may be cumbersome to determine Sp± and A± in (14) 
explicitly. A simpler way is to base the calculations on the 
individual pseudoranges. 3y noting that Pi~t> ~ l^u~^si I 
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(neglecting measurement errors) one may determine Sp± nd ± and 
Aind± that fulfill 

Spindi -A±ndi ^ Pi ~b< Spi nd i +A ±ndi , i = l , ... , XX ( 17 ) 

where 

5 Spxndi is the mean value of the maximum and minimum range to 
satellite i that can be found anywhere within the initial 
uncertainty area and where 

Aindi is half the difference between the maximum and minimum 
range to satellite i that can be found anywhere within the 
10 initial uncertainty area. 

The maximum and minimum ranges to satellite i and hence also 
dpindi and A± ndi can be determined quite easily by applying 
simple geometry. Applying the triangle inequality one may then 
obtain 

15 | {pi-b-Sp ±ndi ) - ( Px-b-Spinai) I < | (pi~b-Sp indi ) / + / ( Pi-b-Sp lndl ) j 

£ Aindi A indl , i=2,...,n (18) 

Reformulating (18) and using the first equality in (13), we get 

Spindi- $Pindl ~ A±ndi ~ A indl < ki*R + V±- p X * < Sp in di~ Spindl+A indi + A± ndl 

(19) 

2 0 Possible values of the reconstructed pseudoranges (6) are 
stepped with a spacing of R meters. The width of the 
uncertainty region is 2 (A± ndi -f-A± ndl ) . So whenever 

2(4ndi +A indl )>R (20) 



there is a risk that the reconstruction cannot be done 
25 unambiguously. 
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A further simplifying expression may be obtained by using only 
the relation (9) and the triangle inequality (in the form //a/- 
/B/ / < jA+Bj) . First note that we may write 

I (Pi-b-pi' - (px-b-p! ' )) l<\Pi-b-p±' l + \pi-b-p x '/ = 
5 / lx u -x S ij-fx u0 -x si l I + 1 1 x u -x sl I - J x u0 -x sl 1 1 < 

jx u -x si -x u0 +x $i l + jx u -x sl -x u0 +x sl I =]x u -x u0 1 + J x u -x u0 1 < 2A (21) 

wherein px ' is the predicted pseudorange to Space Vehicle #i, 
p± ' =lx sl -x u0 l . 

Using (13) and inserting the expression (6) we end with 
10 - 2A < ki+R+Vi-pi'-ip^-px') <2A (22) 

Possible values of the reconstructed pseudoranges (6) are 
stepped with a spacing of R meters. The width of the 
uncertainty region is 4 A So whenever 

4A >R (23) 

15 there is a risk that the reconstruction cannot be done 
unambiguously- Thus, there is a risk for ambiguous pseudorange 
reconstruction when the uncertainty A is more than approxi- 
mately 7 5 km and ranging signal timing measurements are 
peirformed using a truncation interval corresponding to one C/A 

20 code period. 

The problem of ambiguous pseudorange reconstruction due to 
trxincated measurement of ranging signal timing in combination 
with too large uncertainties in apriori location could be 
addressed by having a mobile station measure ranging signal 
25 timing for each ranging signal without any truncation. This 
would however.- require decoding Time Of Week information on each 
measured ranging signal which significantly increases the 
processing delays and may also decrease the detection 
sensitivity of the GPS receiver integrated in the mobile station 
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since it is signif icantly more difficult to decode the Time Of 
Week information than to detect the C/A code boundaries. For 
mobile assisted AGPS, it would also require modifications of the 
messages used to report the ranging signal measurement results 
5 in order to report not only truncated timing but complete timing 
of the ranging signals . 

The present invention addresses the above elaborated problem by 
providing ways of significantly reducing the risk that, in the 
context of AGPS (both mobile station based and mobile station 
10 assisted AGPS) , the apriori location -uncertainty of a mobile 
station causes ambiguous pseudorange reconstruction. At the same 
time the invention also avoids the need for decoding Time Of 
Week on each measured ranging signal and the associated 
disadvantages . 

15 Fig. 5 illustrates a basic method according to the invention for 
determining a location at which ranging signals from at least 
three satellites have been received, wherein said location is 
apriori known to be located within an initial location 
uncertainty area . 

20 The basic method includes performing at least one pseudorange 
selection cycle. All embodiments of methods according to the 
invention includes performing an initial selection cycle 511 
(comprising substeps 501-504 in Fig. 5) . Some embodiments of the 
invention may optionally include one or more additional 

25 selection cycles 512 (comprising substeps 506-508) . 

After performing pseudorange selection in the at least one 
pseudorange selection cycle, said location is determined at step 
505 by utilizing pseudoranges selected after performing said 
pseudorange selection cycles. Thus step 505 may in some 
3 0 embodiments be performed immediately after the initial selection 
cycle 511, and may in other embodiments be performed after one 
or more additional selection cycles 512. 
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The initial selection cycle 511 includes substeps 501-504. 

At substep 501, a first pseudorange Pi* is determined which with 
respect to said initial location uncertainty area is an 
admissible pseudorange associated with a first satellite . 

. 5 At step 502, all pseudoranges associated with at least two 
additional satellites which combined with said first 

pseudorange form admissible relative pseudoranges are 
determined. Mathematically, a pseudoranges p±* associated with a 
satellite i forms an admissible relative pseudorange when 

10 combined with the first pseudorange Px* if expression (14) is 
satisfied. 

At substep 5 03 a set of pseudorange vectors representing all 
possible combinations of the determined pseudoranges associated 
with said first and at least two additional satellites are 
15 formed. Each pseudorange vector thus includes the first 
pseudorange associated with the first satellite (e.g. SV1 in 
Fig. 1) and. a pseudorange associated with each different 
satellite among the at least two additional satellites (e.g. 
SV2 and SV3 In Fig. 1) . 

2 0 At substep 5 04 a set of selected pseudorange vectors is formed 

by selecting at least one vector from said set of pseudorange 
vectors, wherein at least when said set of pseudorange vectors 
includes plur*al vectors, said selecting includes evaluating each 
vector in said set of pseudorange vectors according to a 
25 predetermined rule for initial vector selection. In some 
embodiments of the invention, such evaluation may be performed 
regardless of the number of vectors in the set of pseudorange 
vectors, i.e. also when there is a single vector . In other 
embodiments of the invention, such evaluation may be performed 

3 0 only when there are plural vectors in the 'set of pseudorange 

vectors, i.e. if there is a single vector in the set of 
psedurange vectors, this single vector may be selected without 
performing siach evaluation. The predetermined rule for initial 
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vector selection is selected so as to eliminate pseudorange 
vectors representing less probable combinations of pseudoranges. 
Different embodiments of the invention may apply different rules 
for initial vector selection. For embodiments where the 
5 additional at least two satellites equals exactly two 
satellites, i,e. when each vector . includes pseudoranges from 
three different satellites, the rule for initial vector 
selection would typically be based on comparing the locations 
associated with the different vectors with the initial 

10 uncertainty area. For embodiments where the additional at least ^\ 
two satellites includes more than two satellites, i.e. when each 
vector includes pseudoranges from four or more satellites, the 
alternatives for defining different rules of initial vector ^ 
selection includes also calculating a minimum loss function 

15 value for each new pseudorange value and selecting the 
pseudorange vector having the least minimum loss function value 
or all pseudorange vectors having minimum loss function values 
below a threshold. 

The optional one or more additional selection cycle 512 includes 
20 substeps 506-508. 

At substep 50 6, all pseudoranges associated with at least one 
additional satellite which combined with said first pseudorange ^ 
form admissible relative pseudoranges (i.e. which satisfy 
expression (14)) are determined. 

25 At substep 5 07 a new set of pseudorange vectors are formed 
representing all possible combinations of the set of vectors 
selected in the previous selection cycle and said determined 
pseudoranges associated with said at least one additional 
satellite. Since each pseudorange vector in the new set of 

30 pseudorange vectors include a pseudorange associated with each 
of said at least one additional satellite, the dimension of -the 
new pseudorange vector is increased as compared to the dimension 
of the pseuodrange vectors selected in the previous selection 
cycle . 
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At substep 508, a new set of selected pseudorange vectors is 
formed by selecting pseudorange vectors from the new set of 
pseudorange vectors, wherein at least when said new set of 
pseudorange vectors incluc3.es plural vectors, said selecting 
includes evaluating each vector in said new set of pseudorange 
vectors according to a predetermined rule for subsequent 
selection. In some embodiments of the invention, such evaluation 
may be performed regardless of the number of vectors in the new 
set of pseudorange, i.e. also when there is a single vector. In 
other embodiments of the invention, such evaluation may be 
performed only when there are plural vectors in the set of 
pseudorange vectors, i.e. if there is a single vector in the new 
set of psedurange vectors, this single vector may be selected 
without performing such evaluation. The predetermined rule for 
subsequent vector selection is selected so as to eliminate 
pseudorange vectors representing less probable combinations of 
pseudoranges . Different embodiments of the invention may apply 
different rules for subsequent vector selection, e.g. based on 
comparing the locations associated with the different vectors 
with the initial uncertainty area, based on calculating a 
minimum loss function value for each new pseudorange value and 
selecting the pseudorange vector having the least minimum loss 
function value or all pseudorange vectors having minimum loss 
function values below a threshold. 

Embodiments of the invention could be implemented both in 
control nodes on the network side (e.g. for supporting mobile 
station assisted AGPS) as well as in mobile stations (e.g. for 
supporting mobile station based AGPS) . 

First exemplary embodiments of a method and an apparatus 
according to the invention implemented in the location server 
node 101 of Fig. 1 are illustrated in Figs 8-10. Before 
elaborating on said firrst exemplary embodiments of the 
invention, some details of the structure and processing 
performed in the mobile station MSI of Fig. 1 in connection with 
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mobile assisted AGPS are elaborated in connection with Fig. 6 
and 7 below. 

Fig. 6 is a block diagram illustrating the structure of the 
mobile station MSI. The mobile station MSI includes a cellular 
5 communication module 601, a positioning module 602; a GPS RF 
front end 603, an antenna 604 for communication with the 
cellular network and a GPS antenna 605. The positioning module 
602 includes a CPU 612, memory 610 and a Digital Signal 
Processor (DSP) 611. The cellular communication module 601 

10 wirelessly receives assistance data from the cellular network 
and wirelessly transmits measurement results to the cellular 
network via base stations in the cellular network. The 
assistance data could consist of ephemeris and clock corrections 
for visible satellites, an approximate location of the mobile 

15 station MSI and an approximate GPS system time. Alternatively 
the assistance data could contain explicit assistance data 
intended only for assisting the correlation processing. The 
communication module 601 forwards received assistance data to 
the positioning module 602 using the interface 606 while 

20 measurement results are provided from the positioning module 602 
to the communication module 601 using interface 613 . The 
communication module 601 also provides the GPS RF front end 603 
and the positioning module 602 with a clock reference 607. The 
GPS RF front end module 603 is controlled by the positioning 

25 module 602 using interface 608. 

Fig. 7 illustrates processing performed by the mobile station 
MSI when it receives a positioning request. 

When the positioning module 602 receives a positioning request 
from the communication module 601, it requests the GPS RF front 
3 0 end 603 to provide GPS signal samples at step 7 01. The GPS RF 
front end 603 receives the GPS frequency band through the 
antenna 605, downconverts the signal to baseband, separates the 
signal into in-phase (I) and quadrature (Q) components, samples 
and converts the signals into digital format, and outputs these 
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to the positioning module 602 through interface 609- The 
positioning module 602 stores the received I and Q data in 
memory 610. 

Step 702 define processing performed on each individual ranging 
signal RS1-RS4 which is included in the measurement report 
transmitted at step 707. Please note that even though Fig. 7 
illustrates sequential processing (see step 704) of each 
individual ranging signal, processing related to different 
ranging signals are preferably performed in parallel. 

A ranging signal y as a function of time t received from an 
arbitrary Space Vehicle SV1-SV4 by thie mobile station MSI can in 
a simplified way be written: 

y(t) = a-c(t-x) -d(t-x) *exp{i- (o>ot+o>at+"<J))}+e(t) 

Here a is the amplitude of the received signal, c(t) is the C/A 
code of the Space Vehicle and d(t) is the navigation data bit 
stream (see Fig. 2) . The term x is the unknown delay of the 
signal which is a function of the distance from the Space 
Vehicle to the position of mobile station MSI , co 0 is the GPS 
carrier frequency, co d is the Doppler frequency of the signal, § 
is an unknown phase and e(t) noise. 

At step 7 02 the C/A code boundaries of a ranging signal are 
determined by the Digital Signal Processor 611 in the 
positioning module 602 using correlation that test all possible 
code phase and Doppler shifts for th.e ranging signal. Once step 
702 has been initially completed for a ranging signal, the DSP 
611 maintains synchronization witti said ranging signal by 
tracking changes in C/A code boundary timing of said ranging 
signal . 

If more ranging signals need to be acquired (an alternative YES 
at step 704) step 702 is repeated for* a next ranging signal (as 
already discussed, the processing of step 702 is preferably 
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performed in parallel for several ranging signals and not 
sequentially as indicated by Fig. 7). The decision in step 7 04 
on whether more ranging signals should be acquired or not could 
be based on the number of ranging signals acquired so far (at 
5 least 3 or preferrably 4 ranging signals should be acquired, but 
acquiring more ranging signals would improve the precision of 
the calculated position) and timing requirements (the response 
time for providing a measurement report signal could be 
configured by a parameter to e.g. within 16 seconds of receiving 
10 a positioning request) . ^ 

If enough ranging signals have been acquired (an alternative NO 
at step 7 04) , the GPS Time Of Week (TOW) of a selected point in ^ 
time is estimated at step 705. Please note that preferably, as 
soon as step 702 has been completed for a first ranging signal, 
15 step 705 is performed based on said ranging signal in parallel 
with acquiring additional ranging signals. 

There are several alternatives for how step 7 05 may be 
performed. Typically TOW estimation is based on determining the 
TOW transmitted in the so called Handover Word (HOW) of one 
20 ranging signal (see Fig. 3), preferably the first aquired 
ranging signal, and then compensating for the propagation delay 
from signal transmission by the Space Vehicle until signal 
reception by the mobile station MSI. 

Determining the transmitted TOW can be performed by direct 
25 decoding of the transmitted TOW. This alternative implies that 
data is demodulated at a rate of 2 0 ms and normally requires 
that subframe boundaries are determined followed by decoding of 
the Handover Word, from which the TOW, ie the transmission time 
t ti , can be derived. Each subframe has a length of 6 s, so this 
30 procedure may require that approximately 8 seconds of navigation 
data is collected. TOW demodulation works down to- approximately 
-172dBW, assuming OdB antenna and is in fact the limiting factor 
for detection sensitivity. 
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Alternatively, the transmitted TOW can toe determined by 
reconstruction using correlation techniques. This procedure also 
requires that demodulated data bits are generated, but instead 
of direct decoding, correlation is made with, known transmitted 
5 navigation data bits (e.g. the contents of the so-called 
Telemetry Word and the HOW word which may be sent to the mobile 
station as part of the assistance data) . This requires that the 
GPS time is a priori known to within a few seconds. This 
procedure works to somewhat lower signal levels than direct TOW 
10 decoding, but most likely the performance is limited by the 
tracking loops that may loose lock at such low signal levels. 
Typically phase locked loops or automatic frecjuency control loop 
are employed for this. But it is expected that this will work 
down to say -179dBW. 

15 Compensating for the propagation delay could be performed by 
applying an expected average propagation delay of 77 ms. 
Alternatively a more accurate propagation delay compensation can 
be derived from assistance data received from the cellular 
network by the mobile station MSI according to the principles 

20 elaborated in copending US patent application by inventors Ari 
Kangas and Janos Toth-Egeto filed September 29, 2004. 

At step 706, the positions in time with respect to the selected 
point in time are measured for the next C/A code boundaries 
following the selected point in time. More specifically, for 

25 each acquired ranging signal, the position of the next C/A code 
boundary following the selected point in time is measured by 
registering the number of whole and fractional chips from the 
selected point in time until the next C/A code boundary. 
Finally, at step 707 a measurement report signal reporting the 

3 0 measurement results is wirelessly transmitted by the mobile 
station MSI to the cellular network NETl . The measurement report 
signal is in this example scenario transmitted as ordinary user 
data in the user plane addressed to the location server 101. 
Hence the measurement report signal is transparently routed 
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through the cellular network NET1 via the IP based, network 102 
to the location server 101. 

Fig. 8 schematically illustrates the structure of the location 
server 101 according to a first exemplary embodiment of an 
5 apparatus according to the invention. The location server 
includes a communication module 801 and a positioning module 
803. The communication module 801 receives the measurement 
report and forwards measurement data to the positioning module 
803. The positioning module includes one or more processors CP1 

10 programmed to determine the location of the mobile station MSI 
using the provided measurement data (including the measured 
timing information for each reported ranging signal) and apriori 
information on the mobile station location. The apriori 
information can e.g. be derived from a Public Land Mobile 

15 Network (PLMN) identity included in the signal from the mobile 
station MSI and indicating in which network the mobile station 
MSI is operating. The PLMN identity could e.g. be included as 
part of the cell identity of the cell in which the mobile 
station is currently operating. Using the provided PLMN 

2 0 identity, the positioning module 803 could derive the apriori 

location information e.g. by retrieving the coordinates of the 
centre of the country in which the mobile station MSI is 
operating, together with a radius corresponding to the maximum ^ 
distance from said centre until the border of said country from 
25 a table. The hierarchical nature of cell identities could also 

be exploited, in particular for large countries, to identify a ^ 
particular region within a country in which the mobile station 
is operating. Maintaining a table of centre/ radius information 
for different countries, or regions within said countries, is 

3 0 significantly less burdensome than trying to maintain a global 

data base with information on the geographical coordinates of 
each cell . 

Before discussing more in detail the processing performed by 
- the location server 101 according to an exemplary embodiment of 
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a method according to the invention, some basic calculations 
for determining a mobile station position are described below. 

The measurement equation (2) may be expressed in vectorized 
form as 

5 p =/l n - x u -X s l+bl n +e (24) 

Here p is the column vector of length n containing the 
pseudoranges, l n is a column vector of length n containing only 
ones. X s is a matrix where the ith row contains the coordinates 
x sl of the ith satellite (Space Vehicle) . It is assumed here 
10 that the norm \z\ is calculated for each row of the matrix Z 
within brackets. 

A Taylor series expansion around the inital estimate of the 
unknown parameters x u and b may be expressed as 

p =jl n ' x u -Xsl+(b-b 0 )ln+e =lln m x u0 ~X s j +G • ( (Xu-x u o) (b-b 0 ) ) T +v (25) 

15 wherein G is the geometry matrix which contains the derivatives 
of the pseudoranges with respect to the parameters x u and b, v 
is the sum of the measurement error term e and the higher order 
Taylor series terms. Let 

r± = j x u0 ~x sl j (26) 

20 then G is a matrix with the ith row being equal to 

G± = [ (x uo-x S i) /r± (y u o-ysi)/r± (z u0 -z si ) /r ± 1] (27) 

The least squares solution to (25) is equal to 

( x u b) T =( x u0 b 0 ) T +(G T G)- 1 G T (p-ll n ^ x u o-X s l-boln) (29) 

The least squares minimum loss function value, is equal to . 

25 V=<p-ll n - x u0 -X s l-b o ln) T (T-G(G T G)- 1 G T ) (p-jln' x u0 -X s / -b 0 l n ) (30) 
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The minimum loss function value is a measure of how well the 
predicted pseudoranges based on an vapdated parameter estimate, 
match the measured pseudoranges. 

Fig. 9A-B illustrates processing performed by the location 
5 server 101, and more in particular the positioning module 803, 
according to a first exemplary embodiment of a method according 
to the invention. 

The procedure starts at step 90L by calculating a first 
pseudorange Pi* for an arbitrary first satellite (Space 

10 Vehicle) according to equation (12b) . Further all pseudoranges 
p 2 * (for a second satellite) and p& * (for a third satellite) 
which combined with the first pseudorange pi* would form 
admissible relative pseudoranges are determined using any one 
of expressions (15), (19), or (22). Assuming n 2 and n 3 

15 admissible pseudoranges were found for the second and third 
satellite respectively, then a sest of n 2 * n 3 candidate 
pseudorange vectors can be formed representing all possible 
combinations of the determined pse~udoranges Pi* , p2* and p 3 * 
associated with the first, second and third satellite 

20 respectively. 

At step 9 02 a least squares estimate according to (29) is 
performed for each candidate pseudorange vector, but since only 
three measurements are used, only a two dimensional location 
fix can be made. It is most natural to estimate the x 
25 (latitude) and y (longitude) coordinates since the apriori 
uncertainty in the vertical dimensions is typically smaller than 
the horizontal uncertainty. The reduced parameter vector (jc u y u 
Jb) T is estimated for all candidate pseudorange vectors. 

At step 903 a selection of the most: likely pseudorange vector 
30 candidate (s) is done by considering the size of the update step 
in the horizontal dimension, i.e. Joy comparing the positions 
associated with the candidate pseudorange vectors and the 
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initial location uncertainty area defined by . the apriori 
location x u o of the mobile station and the uncertainty A. Hence 
( (x u - x u o ) 2 +(Yu - Yuo) 2 ) 1 ' 2 is determined for all candidate 
vectors. Since it is known apriori that ( (x u - x u0 ) 2 +(¥u - 
Yuo) 2 ) 1/2 < A for the true parameter vector x u , . y u those 
candidates that generate an update step larger than a defined 
threshold value, selected as larger than A, are excluded from 
further calculations. The threshold could preferably be 
selected larger than A with a margin accounting for the fact 
that the least squares solution normally has to be performed 
iteratively, since the least squares equation (29) is only an 
approximation for the original non-linear problem (24) . 

Note that in this case the residual V in (30) cannot be used, 
since it will in most cases be equal to zero in the case of 
three measurements, three equations . 

In the next step 904 it is checked whether all pseudoranges 
have been reconstructed, i.e. whether all measured satellite 
ranging signals have been considered. 

If all pseudoranges have been reconstructed (an alternative YES 
at step 904) , processing continues at step 913 where the mobile 
station position is calculated/estimated. 

Otherwise (an alternative NO at step 9 04) , processing continues 
at step 905 where all pseudoranges p 4 * associated with an 
additional (fourth) satellite and which combined with the first 
pseudorange Pi* would form admissible relative pseudoranges are 
determined using any one of expressions (15), (19), or (22). A 
new set of all possible candidate pseudorange vectors is formed 
by combining the set of pseudorange vectors selected at step 
903 with all determined pseudoranges p 4 * associated, with the 
additional (fourth) satellite . 
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At step 906, the loss function value (30) is calculated for all 
candidate pseudorange vectors, assuming also in this case that 
only a two -dimensional location vector is estimated. Note that 
in this step no explicit solution ( 29) needs to be calculated. 

5 At step 907 , the candidate pseudorange vector that minimized 
the loss function (30) is selected. 

At the next step 908 it is checked whether all pseudoranges 
have been reconstructed, i.e. wheizher all measured satellite 
ranging signals have been considered. 

10 If all pseudoranges have been recoastmcted (an alternative YES 

at step 908), processing continues at step 913 where the mobile ^ 
station position is calculated. 

Otherwise (an alternative NO at step 908) , in the following 
steps, enough ranging signal measurements are available to 

15 enable a three-dimensional solxition. The procedure for 
pseudorange reconstruction is thus repeated at step 909 for p k * 
where k=5, i.e. all pseudoranges ps* associated with an 
additional (fifth) satellite which combined with the first 
pseudorange p j * would form admissit>le relative pseudoranges are 

20 determined using any one of expressions (15) , (19) , or (22) and ^ 
a new set of all possible candidate pseudorange vectors is 
formed by combining the set of candidate pseudorange vectors 
selected at the previous selection step with all determined ^ 
pseudoranges p 5 * associated wifc-h the additional (fifth) 

25 satellite. 

The loss function values are computed at step 910 for all 
candidate pseudorange vectors in "the new set using (30) . The 
most unlikely candidate pseudorange vectors are excluded from 
further consideration at step 911 based on said loss function 
30 values, i.e. a new set of selected candidate pseudorange 
vectors is formed by excluding unlikely candidate pseudorange 
vectors formed at step 909. In thzls embodiment, the candidate 
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pseudorange vector that minimized the loss function is 
selected. 

At step 912 it is checked whether all pseudoranges have been 
reconstructed, i.e. whether all measured ranging signals have 
5 been considered. If all pseudoranges have been reconstructed 
(an alternative YES at step 912) , processing continues at step 
913 where the mobile station position is calculated. 

Otherwise(an alternative NO at step 912), the procedure of 
pseudorange reconstruction and exclusion of unlikely 
10 pseudorange vectors is repeated from step 909 to determine 
pseudoranges p^* associated with additional satellites k where 
k is incremented, before each new iteration of steps 909-911. 

When all pseduoranges have been reconstructed, i.e. all 
acquired satellite ranging signals have been considered, the 
15 mobile station position x u is calculated/estimated by 
determining the least squares solution (29) at step 913 
possibly using an iterative scheme. 

The pseudorange selection performed in the procedure illustrated 
in Fig. 9 is organized in an initial pseudorange selection cycle 
20 including steps 901-903 and additional pseudorange selection 
cycles 905-907 (operating on pseudoranges associated with four 
satellites) and 908-911 (operating on pseudoranges associated 
with more than four satellites) respectively. 

Fig. 10 illustrates an example scenario of results and 

2 5 intermediate results when performing the method according to 

Fig. 9A-B. In this example scenario, it is assumed that ranging 
signal measurements associated with the four satellites of Fig. 
1 are available, which implies that only the initial selection 
cycle (steps 901-903) and the additional selection cycle (905- 

3 0 907) operating on pseudoranges associated with four satellites 

are performed. 
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Fig. 10 illustrates the first pseudorange p 2 * (rll in Fig. 10) , 
pseudoranges p 2 * (r21, r22 in Fig. 10) and p 3 * (r31, r32 in Fig. 
10) associated with the second satellite SV2 and the third 
satellite SV3 respectively which in this example scenario are 
5 determined at step 901. Based on the determined pseudoranges, a 
set 1001 of pseudorange vectors representing all possible 
combinations of said pseudoranges are also formed at step 901. 

From the set 1001 of pseudoranges formed at step 901 , a set 1002 
of selected pseudorange vectors is derived at steps 902-903 by 
10 selecting pseudorange vectors which are associated with 
locations within the initial uncertainty area. As illustrated in 
Fig. 10, two pseudorange vectors are selected at steps 902-903 _ 
in this particular example scenario. 

Fig. 10 further illustrates the pseudoranges p 4 * (r41, r42 in 
15 Fig. 10) associated with the additional fourth satellite SV4 
which are determined at step 905. Based on the determined 
pseudoranges p 4 * associated with the fourth satellite SV4, a new 
set 1003 of pseudorange vectors representing all possible 
combinations of the pseduoranges p 4 * and the set 1002 of 
20 selected pseudorange vectors from the previous selection cycle 
(i.e. the initial selection cycle in this example scenario) are 
also formed at step 905 . 

From the new set 1003 of pseudoranges formed at step 905, a new 
set 1004 of selected pseudorange vectors is derived at steps 
25 906-907 by selecting the pseudorange vector which minimizes the 
minimum loss function (30) . 

Finally, according to this example scenario, the elements of the 
single pseudorange vector in the new set 1004 of selected 
pseudorange vectors are used at step 913 to determine the 
3 0 location (represented by vector 1005 in Fig 10) at which the • 
ranging signals from the four satellites SV1-SV4 were received. 
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Apart from the exemplary first embodiment of the invention 

disclosed above, there are several ways of providing 

rearrangements, modifications and substitutions of the first 

embodiment resulting in additional embodiments of the invention. 

5 In the first exemplary embodiment of the invention, the 
processing steps illustrated in Fig. 9A-B are performed by 
digital data processing circuitry in the form of one or more 
conventional programmable processors. However, any digital data 

* processing circuitry capable of performing said processing could 
10 be used, e.g. a state machine, an ASIC, a discrete logic circuit 

etc. In the first exemplary embodiment of the invention, as in 

* other embodiments of the invention using programmable devices, 
the controlling computer program (software) is embodied as 
machine-readable instructions stored on some kind of computer- 

15 readable medium such as RAM, a hard drive, electronic read-only 
memory, an optical storage device (e.g. a CD-ROM as 
schematically illustrated in Fig. 11) etc. Programmable devices 
performing processing according to the invention, can be 
dedicated to this task or used also for processing relating to 

20 other tasks. 

An exemplary embodiment of the invention for use in the context 
of mobile station based AGPS could be derived from the 
illustrated first embodiment of the invention by essentially 
replacing step 7 07 of Fig. 7 with processing according to Fig. 

25 9A-9B in the positioning module 602 of the mobilLe station MSI . 
Thus the calculations performed by positioning module 9 03 of the 
location server 101 in the first exemplary embodiment would 
instead be performed by positioning module 602 of the mobile 
station MSI. An apriori estimate of the mobile station position 

30 together with satellite ephemeris data and clock corrections 
would be provided by the network as assistance data for use in 
the position calculations. Thus, the mobile station MSI would 
become an exemplary embodiment of an apparatus according to the 
invention in this scenario. 
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The invention could of course be applied both in the context of 
control plane and user plane solutions to AGPS. As regards 
control plane solutions to AGPS, the invention is probably most 
interesting to apply in the context of extended range cells (in 
5 GSM, extended range cells could have a radius of up to 100 km) 
or when cell identity positioning (which typically is used as a 
basis for determining apriori location information) is not 
implemented in a network. 

An alternative approach to performing pseudorange selection 

10 cycles 909-911 until all measured ranging signals have been 
considered would be to quit the loop of iterative selection 
cycles after a selection cycle resulting in the selection of a 
single pseudorange vector. An updated location uncertainty area 
within the initial location uncertainty area defined by the 

15 apriori location x u0 of the mobile station and the uncertainty 
A could then be determined based on the selected single 
pseudorange vector. The unknonwn location of the mobile station, 
i.e. where the satellite ranging signals were received, could 
then be determined using the updated location uncertainty area. 

20 Since the updated location uncertainty area would be 
significantly smaller than the initial location uncertainty 
area, there is in practice no risk of ambiguous pseudorange 
reconstruction when determining the unknown location using the 
updated location uncertainty area. The actual location 

25 calculations could be performed as for normal GPS-positioning. 

In situations where it would be desirable to handle even larger 
apriori location uncertainties in connection with AGPS, the 
present invention could be combined with the teachings of 
International Patent Application No PCT/IB2004/052040 by 
30 performing measurements modulo the navigation data bit length on 
received ranging signals as specified in said International 
Patent Application application and then determining the location 
at which said ranging signals were received as disclosed in this 
patent application . 
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In order to provide increased robustness against e.g. errors in 
estimated GPS time at signal arrival, position calculations may 
be performed in accordance with copending US patent application 
No 60/545194 by Ari Kangas . 

5 There are several alternative rules that different embodiments 
of the invention may apply when evaluating and selecting 
pseudorange vectors in additional pseudorange selection cycles. 
Instead of applying the minimum loss function (30) , a weighted 
minimum loss function 

10 V = (p-\l n - x u -X 5 j-b 0 l n ) T QT 1 (p-\l n ' Xu -X s l-boln) (31) 

could be applied where Q is the covariance matrix of range 
signal measurement errors . 

Instead of selecting the pseudorange vector which minimized the 
loss function as in steps 907 and 909 in Fig. 9A-9B, all 

15 pseudorange vectors that produce a loss function less than a 
defined threshold could be selected. Further, the evaluation 
rule applied may imply different evaluation criterias in 
different selection cycle iterations. One such example would be 
to select the pseudorange vector which minimized the loss 

20 function in the last iteration, but selecting all pseudorange 
vectors producing a minimum loss function value below a defined 
threshold in all but the last iteration. Another example of 
different evaluation criterias in different cycles would be to 
apply different threshold values for different iterations 

25 (generally applying narrower and narrower threshold values) . 

Evaluation of pseudorange vectors could also be performed by 
calculating locations associated with each pseudorange vector 
and comparing the calculated locations to the initial location 
uncertainty area (i.e. similar to the comparision performed in 

3 0 the initial pseudrange selection cycle) . It is also possible to 
evaluate the pseudorange vectors by a combined criteria based on 
both minimum loss function values (according to e.g equation 
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(30) or (31) ) and a comparision of associated locations and the 
initial uncertainty area. 

In the first exemplary embodiment of the invention, truncated 

timing measurement results are available for all ranging signals 

5 RS1-RS4. In other embodiments of the invention, one or more 

» 

ranging signals could have been measured without truncation. 
However, as long as truncated rangin signal measurements are 
performed for at least one of the satellites, i.e. the available 
measurement results include truncated timing measurement results 
10 for at least one satellite, there is a risk for ambiguous 
pseudorange reconstruction if the apriori location uncertainty 
area is large and hence the present invention may be used to 
provide increased robustness against ambiguous pseudorange 
reconstruction . 

15 Even though the invention in its first exemplary embodiment has 
been applied in the context assisted GPS, the invention may of 
course be applied in connection with other satellite based 
positioning systems such as GALILEO or GLONASS to address 
situations where there is a risk that truncated ranging signal 

2 0 measurements in combination with a large apriori location 
uncertainty area implies a risk for ambiguous pseudorange 
reconstruction . 
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